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. ' Abstract. We perform Ramsey-Borde spectroscopy on laser cooled magnesium 

fj . atoms in free fall to measure the ^^o — > '^Pi intercombination transition frequency. 

The measured value of 655 659 923 839 730 (48) Hz is consistent with our former 
^~»' atomic beam measurement [1]. We improve upon the fractional accuracy of the 

i-C , previous measurement by more than an order of magnitude to 7 x 10~^^. The 

magnesium frequency standard was referenced to a fountain clock of the Physikalisch 
Technische Bundesanstalt (PTB) via a phase-stabilized telecom fiber link and its 
stability characterized for interrogation times up to 8000 s. The high temperature 
of the atomic ensemble leads to a systematic shift due to the motion of atoms across 
IQ . the spectroscopy beams. In our regime, this leads to a counterintuitive reduction of 

residual Doppler shift with increasing resolution. Our theoretical model of the atom- 
light interaction is in agreement with the observed effect and allows us to quantify its 
contribution in the uncertainty budget. 
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Remote frequency measurement of the ^Sq — )■ ^Pi transition in laser cooled "^^Mg 2 

1. Introduction 

The inetastable ^Pi state in magnesium has a lifetime of 4.4 ms [2] which is longer than 
that for heavier alkaline earth elements. This enables us to perform high resolution 
Ramsey-Borde interferometry [3], a well established technique in spectroscopy [4, 5, 6, 7], 
for interrogation. We apply this method to determine the frequencies of optical 
transitions from the ground state to the metastable states ^Pi and ^Pq. The frequency 
value of the latter transition is derived by combining our measurement for the ^Pi state 
with measurements from [8]. This transition is of relevance for the implementation of 
a lattice clock [9], which is currently underway in our laboratory. Magnesium offers 
an advantage in that the transition is much less perturbed by blackbody radiation 
compared to heavier alkaline earth elements currently used for optical clocks [10, 11]. In 
strontium, for example, one way to minimize this effect is by performing measurements 
in a cryogenic environment [12]. 

The frequency measurement is carried out on an atomic ensemble at a temperature 
of few milliKelvin and with size comparable to that of the interrogation beams. 
This makes the measurement very sensitive to phase shifts induced by the wavefront 
curvature, which are also of concern in atom interferometers [13, 14]. These effects 
can be minimized by lowering the temperatures, thus yielding more compact clouds, 
and larger interrogation beam diameters, neither of which were feasible in our case due 
to high Doppler temperature and laser power constrains. Atom interferometers employ 
methods such as using interrogation beams in a retro- reflecting configuration [15], which 
cannot be applied to a Ramsey-Borde interferometer used in this work. 

The first order Doppler effect caused by the milliKelvin atomic ensemble in 
combination with the finite wavefront curvature led to a different behavior in our 
measurement than previously reported by Wilpers, et al. [16] in the context of the 
Calcium atomic clock. In our case, the atoms explore a large fraction of the interrogation 
beam in the transverse direction. This results in atoms on the edge contributing 
differently, both in terms of phase and excitation probability, to the Ramsey interference 
pattern than those at the center. We extend the work in [16] to develop a theoretical 
model to account for the effects that occur in our regime. The model is described in 
section 3.1 in the context of the uncertainty budget of the frequency measurement. 

2. Experimental apparatus 

2.1. The Ramsey-Borde interferometer 

The frequency measurements were performed by locking an ultrastable laser to a time- 
domain Ramsey-Borde atom interferometer [3, 4]. Figure 1(a) shows the partial ^^Mg 
level scheme with all relevant transitions. Magnesium is heated in an oven to 420° C to 
produce an effusive thermal beam. Laser cooling is carried out on the ^5*0 — )■ ^Pi closed 
transition at 285 nm and the atoms are captured in a magneto optical trap (MOT). The 
light for driving this transition is produced by a Rhodamine 6G dye laser which outputs 
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Figure 1. (a) The electronic states of ^^Mg and tlie wavelengths of relevant 
transitions used to manipulate atoms, (b) Temporal sequence of the experiment 
includes preparation, coherent manipulation (spectroscopy) and detection of atoms. 
The spectroscopy part of the sequence is shown in detail. Arrows indicate the direction 
of laser beams. Tq (= 1-5 ms) is the delay at the start of the interferometer after 
preparation, Tp is the duration of the light pulses and T' (= 20 /is) is the time between 
the two pairs of light pulses. The coherent evolution time between the first and the 
second light pulse is T. The four pulse Ramsey-Borde interferometer sequence itself 
lasts for 2T + 4:Tp+T' sa 0.5 ms. The interferometer is always operated in a symmetric 
configuration where the delay between the third and fourth light pulse is also T. For 
the frequency measurement, T = 204 /xs and Tp = 3.0 /iS. 



1.2 W of power at 571 nm, followed by cavity enhanced second harmonic generation 
(SHG) inside a /3-BaB204(BBO) crystal which yields 60 mW at 285 nm. Light from the 
dye laser is transported to the SHG via a polarization maintaining fiber with ^ 50% 
efficiency. The broad linewidth (78 MHz) of the cooling transition implies a relatively 
short Zeeman slower length of the order of 10 cm. This small distance allows us to use 
the quadrupole magnetic field of the MOT itself (0.65 T/m along the radial direction) 
as field for the slowing, avoiding additional coils. The MOT field together with a laser 
beam carrying 15 mW of power, red detuned by 2.5 full linewidths from resonance and 
counter propagating to the atomic beam increases the captured atom number by a factor 
of 25. We load 2.5 x 10^ atoms in the MOT in 250 ms at a peak density of 4 x 10^^ m~^. 
The ensemble has a temperature of 3 mK (which is slightly above the Doppler limit of 
1.9 mK) for a laser detuning of one full linewidth below resonance. 

After the loading sequence, the cooling light and the MOT magnetic field are turned 
off within 700 /is. Atoms are released from the trap at this stage and the cloud expands 
ballistically in a homogeneous magnetic field of 0.6 mT. This magnetic field is created 
by large Helmholtz coils and is essential to define the quantization axis allowing the 
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Figure 2. (Left) Typical interference scans obtained by tuning the laser frequency. 
The frequency detuning is with respect to the center of the atomic transition. Data 
points arc individual runs without averaging. The fringe period and hence the 
frequency resolution was adjusted by changing the time T. T (= 40.6, 128.6, 204, 405 
[is) increases from top to bottom. The pulse duration Tp = 3.0 \xs and time between the 
laser pulse pairs T' = 20 [is. (Right) Fractional instability of the frequency standard 
as a function of resolution obtained from the Ramsey pattern according to equation 1. 
The solid line is the theoretical stability for a signal to noise ratio of 3. 



separation of different Zeeman sublevels of the ^P\ state. 

The interferometer is operated on the ^5*0 — ?■ ^P\ transition and the mj = — ?■ mji = 
clock transition is probed. The temporal sequence of the experiment is depicted in 
figure 1(b). The interferometer sequence starts 1.5 ms after the release when effects due 
to the decaying magnetic field are negligible. The peak density drops to 2 x 10^^ m~^ 
due to the expansion. Four 7r/2 pulses of duration Tp are used to coherently split and 
recombine the atomic ensemble. The direction of the second pulse pair is reversed with 
respect to the first one. The time T' between first and second pair of pulses is fixed 
to 20 yus. The experimental sequence is dominated by the time for preparation and 
detection of atoms. The evolution of atoms in the Ramsey-Borde interferometer lasts 
for 2T + 4Tp+T' ^0.5 ms. The time T between the first and the second pulse and 
Tp both determine the frequency resolution A of the interferometer and is given by 
A = 1/(4T + 4:Tp) [17]. The spectroscopy (clock) laser is an amplified diode laser at 
914 nm that is stabilized to an ultrastable high finesse resonator [18]. A fiber frequency 
comb is used to compare the frequency of the clock laser to a reference frequency source. 
The stabilized clock laser output is frequency doubled using cavity enhanced SHG in 
a non-linear KNbOs crystal and used for interrogation. The interrogation beams are 
coUimated to a 1/e^ radius of 2 mm. 

The large Doppler width of the ensemble (3 MHz), small homogeneous width of 
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the 7r/2 pulse (0.3 MHz) and limited power in the interrogation laser (40 mW) limits 
the fraction of excited atoms to 1%. In order to detect this small fraction, we measure 
the excited state population by pumping those atoms to the ^P2 state via the ^D2 state. 
The atoms are then captured in a MOT based on the ^P2 — ?■ ^D^ cycling transition 
in the metastable manifold. The relatively smaller fine structure splitting between the 
different ^D states results in off resonant excitation to ^Di and ^Z)2 states. These states, 
after spontaneous decay, populate the ^Pq and ^Pi states. Hence, detection is done in 
the presence of the repump laser for the ^Pi state and an additional repumper for the 
■^Pq state via the ^Di state. The light for these three transitions is produced by amplified 
diode lasers at 767 nm followed by cavity enhanced SHG in a LiBaOs (LBO) crystal. 
The cooling laser is capable of delivering a total of 15 mW at 383 nm to the atoms 
after losses in the acousto optic modulator (AOM) and fiber coupling. The repump 
lasers can deliver a total of 2 mW each. With this background free detection scheme, 
we greatly increase the signal-to-noise ratio in comparison to conventional ground state 
detection. Flourescence light of the metastable MOT is collected by a lens with a 
collection efficiency of 10~^ and imaged on a photomultiplier tube with a quantum 
efficiency of 25 %. In order to avoid limitations by photon shot noise, we choose the 
detection time to be 50 ms giving a typical photon count of 8 per atom. The cycle time 
of the frequency measurement is 400 ms. 

Figure 2 shows Ramsey fringes obtained for typical experimental parameters. All 
data points were obtained without any averaging. The observed fringes were fitted with 
sinusoidal functions. The amplitude of the fit gives the signal while the root mean 
square of residuals are used to determine the noise. The fractional frequency instability 
at 1 second is estimated using [19], 

1 

7rQ~S/N 



^.(^) = _ ^ o /.T \lTc/r. (1) 



Here, r is the integration time, Q = vq/ /S., vq is the transition frequency, 2A is the 
fringe width, T^ is the cycle time and S/N is the signal-to-noise (i.e. fringe amplitude 
(due to contrast) to rms noise) ratio. Figure 2 also shows the calculated instability 
at 1 s against the interferometer resolution. The highest stability can be expected 
for resolutions finer than 1 kHz. For coarser resolutions, the signal-to-noise ratio is 
constant and the instability increases linearly. In this regime, the stability is limited by 
fluctuations in the atom number. For finer resolutions below 1 kHz, the fringe contrast 
drops considerably as is shown in figure 2. We expect vibrations to be the main cause 
for this effect as degradation due to wavefront curvature (discussed in the next section) 
alone is not sufficient to explain this behavior. 

The interferometer signal is composed of two Ramsey fringe patterns separated by 
twice the recoil shift of 39.8 kHz. To obtain the highest contrast, the dark time T is 
ideally chosen in a way to maximize the overlap of the two Ramsey fringe patterns . This 
happens at a discrete set of values for T. The measurements in figure 2 were performed 
only at these specific values. From these measurements we expect the highest stability 
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Figure 3. Fractional Allan deviation of the frequency of the magnesium standard 
compared to different references. The black squares show the instability with respect 
to a local passive hydrogen maser. The red circles (blue triangles) show a remote 
comparison with respect to an active hydrogen maser (ultrastable laser) at PTB. The 
dashed black (red) line is the typical stability of the passive (active) maser. The best 
characterization is given by the open blue circles that corrects for the drift in the 
ultrastable laser for long timescales using the active maser. The measurement time 
covered more than 8000 s allowing us to calculate the Allan deviation for averaging 
times up to 2000 s. During unlocking events of the dye laser, which in total amounted 
to less than 250 s, a linear drift correction was applied to the clock laser. 



for the frequency measurement at T = 204.0 /is which corresponding to a resolution 
of 1.2 kHz. In order to generate an error signal which is insensitive to linear and 
quadratic background offset of the signal, we used a square wave modulation technique 
where the excitation probability is probed at four detunings on the steepest slope of the 
interferogram [20]. 



2.2. Remote stability characterization of the frequency standard 

The magnesium frequency standard was characterized with respect to a passive hydrogen 
maser (H7), stationed locally, and an active hydrogen maser (H5) and the cesium 
fountain clock CSFl [21, 22], located at PTB, via a phase-stabilized fiber link. Details 
of the link can be found in [18]. The transfer of stability via a frequency comb from the 
ultrastable laser to a 1.5 /im laser can be achieved with an instability of 10~^^ in 1 s and 
5 X 10~^^ in 10,000 s [23]. The link does not limit the determination of uncertainty of 
the frequency measurement on any timescale. Even for very long timescales, it does not 
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Figure 4. The previous [1] and current frequency measurement of the ^Sq — > ^Pi 
intercombination transition in magnesium. The 2007 value was obtained with a thermal 
atomic beam apparatus using a transportable Cs atomic clock of the PTB as reference. 
The current value was obtained by performing spectroscopy on laser cooled atoms with 
a Cs fountain clock as reference that was stationed remotely at the PTB. 



introduce any uncertainty above the statistical limit [24]. We estimate the fractional 
uncertainty due to the fiber link to be < 1 x 10"^^ for an averaging time of 1000 s. 

The results of the frequency comparisons are shown in figure 3. The local 
comparison of the clock laser with H7 is dominated by the instability of the microwave 
system which is 4 x 10~^'^/-i/r/s. The comparison with the ultrastable laser at PTB, 
with a typical stability of few times 10~^^, resolves the short term stability of the 
magnesium frequency standard. The slight increase in instability above 10 s is due 
to the tight locking of the clock laser to the atomic resonance. The tight lock allows 
faster convergence of the laser frequency. Beyond 100 s the instability measurement is 
dominated by the drifting ultrastable laser at PTB which is not stabilized to an absolute 
reference. The instability of the magnesium frequency standard for these timescales is 
resolved with H5. For averaging times between 100 s and 1000 s, we observe the stability 
of the magnesium frequency standard to be 2 x 10~^^/wr/s. We can study systematic 
shifts below 10~^^ after an integration time of 500 s . 



3. Frequency measurement 



Several frequency measurements were performed with our frequency standard linked to 
the fountain CSFl at PTB. The measurements result in an intercombination transition 
frequency of 

v{^^M^}Sq -> ^Pi) = 655 659 923 839 730 (48) Hz. (2) 



Remote frequency measurement of the ^Sq — )■ ^Pi transition in laser cooled "^^Mg 8 

Table 1. Uncertainty budget for the frequency measurement on laser cooled 
magnesium atoms. 



Effect 


Frequency correction 


Uncertainty 


Residual first order Doppler 


OHz 


42 Hz 


AOM phase excursion 


OHz 


20 Hz 


Blackbody radiation 


-0.23 Hz 


+0.04 TT„ 

-7.0 -n-^ 


Quadratic Zeeman 


+61 Hz 


7 Hz 


Subtotal 


+61 Hz 


47.6 Hz 


Statistics of frequency measurement 


OHz 


8 Hz 


Gravitational red shift 


+5 Hz 


0.4 Hz 


Fiber link 


OHz 


6 mHz 


Other (from sec 3.2.4) 


OHz 


2.4 Hz 


Total 


+66 Hz 


48.3 Hz 



Figure 4 shows the new frequency value, which is consistent with our previous 
measurement [1]. The applied corrections due to various effects and their uncertainty 
contributions are summarized in Table 1. The total uncertainty of 48 Hz amounts 
to a fractional uncertainty of 7 x 10~^^. The individual effects and their uncertainty 
contributions are discussed in the following sections. 

3.1. Residual Doppler effect 

The relatively high temperature of the freely expanding atomic ensemble results in a 
significant motion of atoms during the interferometer sequence. This motion causes 
a residual Doppler shift and manifests itself via two main effects. Firstly, it causes 
dephasing within the atomic ensemble due to the finite wavefront curvature. Secondly, 
it modifies the excitation probability of atoms within the ensemble due to the Gaussian 
intensity profile of the laser beam. Both these effects are exacerbated by the large initial 
size of the atomic ensemble at the beginning of the interferometry sequence which is of 
the same order as the interrogation beam diameter. 

During each interferometer pulse, an atom samples the local phase and intensity 
of the light field. Wavefront curvature, beam tilts or an inhomogeneous intensity 
distribution combined with the movement of the atom between pulses leads to a 
phase shift in the interferometer. This phase shift enters the interferometric phase 
via A0 = 04 — 03 + 02 — 01! where, 0j denotes the phase of the electromagnetic wave 
at the position of the atom during the i^^ interrogation pulse. For phase shifts A0 ^ vr 
this leads to a frequency shift Az/ = A0/(47r(T + r^)) during the measurement [16]. 

We performed a simulation to calculate the frequency shift due to these effects. 
The simulation takes into account temperature and size of the atomic cloud, intensity 
and diameter of the laser beam, wavefront curvature and beam tilts and displacement 
of the cloud center with respect to the laser beam. The trajectory followed by an 
individual atom results in a Ramsey fringe with amplitude and phase shift given by 
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Figure 5. Schematic illustrating the position dependent phase and intensity of the 
light field experienced by an atom due to its transverse motion during the Ramsey- 
Borde sequence. The interrogation beam wavefront has a radius of curvature R and 
a Gaussian intensity profile. Filled (open) circle represents the position of the atom 
during the first (second) pulse. The wavefront curvature causes a phase shift At/) 
between the two pulses while the Gaussian intensity profile results in a different Rabi 
angle at the end of each pulse. 



the phase and laser intensity at the interrogation pulses, which can be described by a 
complex amplitude. This complex amplitude when integrated over the entire ensemble 
results in the overall frequency shift and fringe amplitude. In addition to the analytical 
treatment of [16], our model considers the different Rabi angles for atoms and is valid 
for large phase shits. These additional effects are important in our regime and result 
in a better agreement with the experiment compared to the analytical treatment alone. 
In the following subsections we describe these effects and compare our simulation with 
experimental data. 



3.1.1. Dephasing Atoms moving away from the center of a curved wavefront see the 
phase growing quadratically with distance (figure 5). For low temperatures the phase 
shift remains below tt. Averaged over the entire ensemble this results in a linear increase 
in frequency shift with respect to T as predicted by the formalism in [16] (dashed lines 
in figure 6(a)). For higher temperatures, the atoms travel further away from the center 
and the interferometric phase begins to lie outside the interval [-7r,+7r]. Hence, the 
Ramsey fringe patterns of atoms with different velocities begin to dephase when the 
complex fringe amplitudes are averaged over the atomic ensemble. As fringe patterns 
from fast moving atoms, that show the largest frequency shift, start to dephase first, 
they contribute less to the overall interference pattern which is now dominated by the 
low velocity atoms that remain close to the center of the cloud. Thus the resulting 
phase shift A0 is bound and the corresponding frequency shift Az/ = A(f)/{An{T + r^)) 
is reduced with increased T (i.e. for finer resolution). In addition, as T increases the 
contrast keeps degrading until the interference pattern is completely washed out. 

A numerical simulation of the frequency shift due to dephasing with respect to 
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Figure 6. (a) Calculated frequency shift versus the time delay T in the Ramsey- 
Bordc sequence. The phase accumulated by atoms in the light field depends on their 
position and velocity. Averaging the phase shifted Ramsey fringes of the individual 
atomic trajectories leads to a signal with a frequency shift relative to the signal from 
atoms at rest. In order to illustrate the different effects, we assume a homogeneous 
intensity profile optimized for a 7r/2 pulse but a curved wavefront with a radius of 
50 m. Dashed lines are calculations according to the model in [16]. (b) Frequency 
shift vs pulse duration Tp for different wavefront curvatures. The temperature of the 
atoms is 3 mK. Here, we assume an inhomogeneous Gaussian intensity profile of radius 
wq = 2 mm. In both figures, the root-mean-square radius of the ensemble upon release 
is set to 750 /im and Tq = 1-5 nis according to experimental values. The calculated 
dots are connected to serve as guide to the eyes. 



pulse separation time T for laser cooled ensembles of magnesium atoms of different 
temperature is shown in figure 6(a). Temperatures of 3 mK, as typically achieved in the 
experiment, along with moderate curvatures (R = 50 m) lead to significant frequency 
shifts of 140 Hz. Only at very short time scales one can treat our scenario with the model 
in [16]. This frequency shift starts to decrease at 400 /is with increasing evolution time T 
and hence with resolution. In order to separate the additional effect due to the Gaussian 
intensity profile of the beams, we assume for this simulation, that the beam diameter 
is much larger than the atomic ensemble. The radius of curvature of the wavefronts is 
R = 50 m and converges along the propagation direction. The two light beams pulsed 
from opposite directions are assumed to have a similar curvature with the same sign 
due to our alignment procedure. The effects due to wavefront curvature of the two light 
fields do not compensate each other in this case. The aforementioned effect will be 
modified by the Gaussian intensity profile of the laser beam, which will be discussed in 
the next subsection. 
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Figure 7. (a) Measured frequency versus the pulse separation T between the first 
and second laser pulse of the Ramsey-Borde sequence for two different pulse durations 
Tp = 3.0, 4.6 us (red squares and blue circles, respectively), (b) Measured frequency 
versus pulse duration Tp for curved wavcfronts with i? = 18 m. Both: The red and 
blue curves show the numerically calculated frequency shifts considering the Gaussian 
shaped laser beam. In order to compare the functional dependence of our model with 
experimental values, we allowed for a frequency offset to account for systematic shifts 
such as a mismatch between the atomic cloud center and the optical axis of the laser 
beams which may cause additional pulse time dependent phase shifts. The red and 
blue horizontal lines show the unshifted frequencies inferred from our model. The 
solid gray lines show the mean value of the measurements with optimized wavcfronts 
and the dashed gray lines indicate their uncertainty as derived from the Monte Carlo 
simulations. The blue data points in (b) have larger error bars because they were 
measured only with respect to a local passive maser. 



3.1.2. Effect of the Gaussian intensity profile The atoms can spread over a large 
fraction of the beam diameter, perpendicular to the propagation direction. Atoms on the 
outer edge of the ensemble are in a low intensity region due to the Gaussian intensity 
profile of the beam (figure 5). This reduces the Rabi angle of these atoms for each 
pulse and they contribute less to the interference pattern. The power of the beams is 
unchanged for all pulse times. For longer pulses, atoms in lower intensity regions in the 
wings have a higher excitation probability and hence contribute more to the ensemble 
average. The frequency shift thus increases. Figure 6(b) shows a simulation of the 
frequency shift with respect to pulse duration. The power is held constant for different 
pulse times which results in a different Rabi angle for the different pulse durations. The 
frequency shift is higher for stronger curvature of the wavefront since atoms experience 
larger phase changes. 



3.1.3. Comparison of simulations with experiments We analyzed the influence of both 
effects by measuring the phase shift as a function of the resolution and the length of 
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Figure 8. (a) Measured frequency values of the ^5*0 -^ ^Pi transition obtained from 
two different sets of wavefront curvatures at a resolution of 1.2 kHz. The error bars only 
account for the statistical errors. All measurements before October 2009 were obtained 
with curved wavefronts and were corrected according to the theoretical model. Blue 
and red data points show measurements with a pulse duration of Tj, = 3.0 jjls and 
Tp = 4.6 /iS, respectively. Black squares are measurements with optimized wavefronts. 
For the optimized setting a pulse duration of 3.0 ^s was applied corresponding to a Rabi 
angle of 7r/2 for the atoms in the center of the beam. The solid gray line is the weighted 
average and the dotted line is the statistical uncertainty, (b) Measured frequency shift 
for two pulse durations Tp = 3.0 /iS and Tp = 4.6 jjls : {^{Tp = 3.0 fj,s) — v{Tp = 4.6 /J.s)). 
Red data points are for curved wavefronts and black data points are for optimized 
wavefronts. The solid red and gray lines indicate the weighted average and dashed 
lines the error bars for the respective data sets. 



the excitation pulses for a wave front curvature of R = 18 m and compared this with 
our model. As the two oppositely directed beams of the Ramsey-Borde interferometer 
were aligned in the same way, they had a similar curvature. This value was determined 
by measuring the beam radius over a long distance. 

Figure 7 shows the results of the first measurement for this relatively strong 
wavefront curvature. Figure 7(a) depicts the variation of the measured frequency with 
respect to the pulse separation T (and hence resolution) for two different pulse durations 
Tp = 3.0, 4.6 /is. The measurements confirm that the frequency shift in this regime 
decreases with increasing resolution. According to our model, the interference pattern's 
offset phase, which shifts, is very sensitive to the pulse duration. Figure 7(b) presents 
the frequency values of the intercombination transition obtained for different durations 
of the excitation pulses. The different interference scans were performed with a constant 
intensity. The scans with different pulse duration were performed in random order to 
minimize other effects such as frequency drifts. Pulse duration of 3.0 fis corresponds 
to a Rabi angle of tt/2 for atoms at the intensity maximum and the Rabi angle scales 
linearly with pulse duration. 

The model correctly predicts the behavior of frequency versus pulse duration for 
data taken on a single day which gives an increased confidence in our model. However, 
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the model accepts many parameters that can fluctuate between measurement runs on 
different days. For example, the frequency values depend critically on the position 
of the atomic ensemble with respect to the spectroscopy beams for the mentioned 
wave front curvatures. The frequency can deviate by as much as 5 Hz for a 1 mm 
misalignment . Based on our model, we can analyze these fluctuations with a Monte 
Carlo simulation for different experimental configurations to determine an uncertainty of 
the absolute frequency shift. In the Monte Carlo simulation the parameter values such as 
temperature, position of the ensemble, beam tilts, etc are varied within their estimated 
uncertainties. This changed the absolute values only by a few percent compared to 
simulations done with only the expectation values of the parameters. We obtained a 
correction of +102 Hz and +67 Hz for frequency values obtained with pulse durations 
of 4.6 /is and 3.0 /xs respectively. 

Figure 8(a) shows the measured frequency of the intercombination line at a 
resolution of 1.2 kHz. A correction has been applied to all data points measured with the 
wavefront of radius R = 18 m. The corrected and uncorrected data points are in good 
agreement, being a further confirmation of our model. The corrected values for curvature 
of radius R = 18 m have large error bars and they make only a minor contribution to the 
final frequency value. Thus, for the determination of the frequency value of the ^5*0 — !■ ^Pi 
transition we only considered the value obtained for R > 30 m. For the optimized wave 
fronts we still observed a frequency shift for different pulse durations, but it was reduced 
by a factor of two on average as can be seen in figure 8(b). We would like to emphasize 
that our model allowed us to calculate accurately the uncertainty contribution due to 
the residual Doppler effect. In order to estimate the uncertainty we assume a wavefront 
curvature of R = 30 m with an identical sign, which is a conservative lower limit for 
the magnitude of the curvature. The uncertainty is obtained by performing a Monte 
Carlo simulation in which we also vary the other experimental parameters within their 
uncertainty. This result gives a mean frequency shift of 41 Hz with a standard deviation 
of 7 Hz. Measurement of such a large radius of curvature gave an uncertainty not only in 
the magnitude but also in the sign of the curvature. This would result in the frequency 
shift predicted by our model to have the same magnitude but opposite sign. Thus, the 
correction to the frequency predicted by our model has to be treated not as a frequency 
shift but as a total uncertainty equal to 42 Hz. 

3.2. Further contributions to the uncertainty budget 

3.2.1. AOM phase excursions A spurious phase shift due to transient ringing in the 
excitation pulse results in a frequency shift as demonstrated in [25]. The bandpass of the 
combined system of switches, RF-amplifiers and the AOM were identified as the reasons 
for this effect. This shift scales inversely with respect to the pulse duration Tp and the 
pulse separation T. From the measurements similar to those shown in figure 7(a), we 
estimate an uncertainty of 20 Hz in the frequency shift. 
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3.2.2. Blackbody radiation The blackbody radiation (BBR) due to finite temperature 
of the surrounding environment shifts the frequency of the intercombination hne. A 
major part of the BBR spectrum does not extend into the optical domain which 
allows one to treat the field in a quasi-static approximation. The mean square of the 
electric field due to BBR at temperature ^ is {E"^) = 8.55 x 10"^ VVlm^K^) x 'd^. 
Since the ^Pi state is anisotropic, the DC polarizability a depends on the direction 
of the applied electric field. BBR contributes equally from all directions and one has 
to calculate the weighted average of the differential polarizability which is given by 
Aa = |a(='Pi,m=±i) + |a(^Pi,m=o)-«(%) = 7.3(0.3) kHz/(kV/cm)' [26, 27, 28, 29, 30]. 
We estimate the temperature of the science chamber to be 293(10) K and that of the 
oven to be 693(10) K which contributes from a solid angle of only 2 x 10"'^. Thus, 
the necessary correction to the frequency is —0.23(0.04) Hz. However, since we do not 
know the exact reflectivity of the chamber, we make a worst case approximation and 
assume that the blackbody radiation from the oven contributes from all directions giving 
a maximum shift of —7.2 Hz. This results in a shift with an asymmetric uncertainty 
contribution of —0.23 (l?;!)^) Hz. Even for the worst case scenario, the effect is of minor 
relevance for the uncertainty budget. 

3.2.3. Zeeman effect The bias magnetic field applied during spectroscopy along with 
the earth's magnetic field and the residual quadrupole field from the MOT cause a shift 
in the clock transition frequency. The largest contribution to this shift comes from the 
second order Zeeman effect. The coupling coefficient for this interaction is 164 Hz/(mT)^ 
[1] . We determined the value of the magnetic field by measuring the frequency shift of 
the mj = ±1 sublevels and obtain values between 11.5 MHz and 14.0 MHz. Using the 
known splitting of ±21 MHz/mT for the mj = ±1 states, we calculate the field to be 
0.61(12) mT. This corresponds to a frequency shift of 61(7) Hz. 

3.2.4. Other contributions Further contribution to the uncertainty budget results from 
a statistical error of the frequency measurement itself and amounts to 8 Hz. Other effects 
that make minor contributions are stated below. The servo error due to imperfectly 
compensated resonator drift causes an uncertainty of 2 Hz. CSFl has an instability 
of 1.4 X 10~^^ at 1 s and a systematic uncertainty of 7.6 x 10"^^ which corresponds 
to 0.5 Hz in the optical domain. The statistical uncertainty in comparing the active 
maser (H5) to CSFl is 1.2 Hz in the optical domain for a measurement time of 5000 s. 
The total uncertainty due to the frequency reference is thus 1.3 Hz. The AC Stark 
shift due to the cooling and spectroscopy lasers causes an uncertainty of 0.4 Hz. The 
gravitational redshift due to difference in height of the magnesium frequency standard 
and the frequency reference at the PTB is +5(0.4) Hz. Second order Doppler effect 
causes a frequency shift of -13(2) mHz. The fiber link causes an uncertainty of 6 mHz. 
DC Stark effect, cold collisions and asymmetric background cause insignificant shifts. 
The total contribution from these effects keeps the frequency value unchanged while 
increasing the uncertainty value by 0.2 %. 
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Table 2. Frequency of the strongly forbidden ^5o -^^Po transition for different isotopes 
of magnesium. The values are calculated using results from this work together with 
data from [8, 31, 32]. 





Isotope 


Frequency (Hz) 


Uncertainty (Hz) 




24Mg 
25Mg 

26Mg 


655 058 646 681 861 
655 060 050 580 000 
655 061 328 149 000 


48 

102 000 

20 000 


4. Conclusion 









We presented the first frequency measurement of the ^5*0 — > ^Pi transition with laser 
cooled magnesium atoms. The frequency measurement and the characterization was 
performed by means of PTB's ultrastable and high accuracy frequency references via an 
optical fiber link. Benefiting from a reduced Doppler shift, we improved the accuracy of 
the transition frequency by a factor of 30 compared to our previous measurement. As 
shown in figure 4, the new value is consistent with our previous measurement. Using 
our accurate spectroscopy measurement and other data in hterature, we calculate the 
frequency of the strongly forbidden ^5*0 — >■ ^Pq transition for all magnesium isotopes 
(Table 2). These transitions are of interest for future optical lattice clocks, one of which 
is under construction in our group. 
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